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Which radio galaxies can make the highest-energy cosmic rays? 
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ABSTRACT 

Numerous authors have suggested that the ultra-high energy cosmic rays (UHECR) detected 
by the Pierre Auger Observatory and other cosmic-ray telescopes may be accelerated in the 
nuclei, jets or lobes of radio galaxies. Here I focus on stochastic acceleration in the lobes. I 
show that the requirement that they accelerate protons to the highest observed energies places 
constraints on the observable properties of radio lobes that are satisfied by a relatively small 
number of objects within the Greisen-Zat'sepin-Kuzmin (GZK) cutoff; if UHECR are pro- 
tons and are accelerated within radio lobes, their sources are probably already known and 
catalogued radio galaxies. I show that lobe acceleration also implies a (charge-dependent) up- 
per energy limit on the UHECR that can be produced in this way; if lobes are the dominant 
accelerators in the local universe and if UHECR are predominantly protons, we are unlikely to 
see cosmic rays much higher in energy than those we have akeady observed. I comment on the 
viability of the stochastic acceleration mechanism and the likely composition of cosmic rays 
accelerated in this way, based on our current understanding of the contents of the large-scale 
lobes of radio galaxies, and finally discuss the impUcations of stochastic lobe acceleration for 
the future of cosmic ray astronomy. 
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1 INTRODUCTION 

It has been known for many years (e.g. Hillas 1984) that the large- 
scale structures of radio-loud active galaxies are possible sites for 
the acceleration of the highest-energy cosmic rays yet to be de- 
tected, the ultra-high-energy cosmic rays (UHECR) with energies 
above a few x 10^® eV. Radio galaxy jets, hotspots and lobes are 
particularly interesting to modellers, both because the synchrotron 
emission by which we see them in the radio already implies the 
presence of a high-energy particle population (albeit leptonic and 
of much lower energies) and therefore of a particle acceleration 
process, and because the physical conditions, in particular the mag- 
netic field strength B, can either be estimated from equipartition 
or minimum energy arguments (Burbidge 1956) or, more recently, 
determined directly from observations of inverse-Compton emis- 
sion (e.g. Hardcastle et al. 2002). It is thus reasonably easy to say 
whether any given component of a radio galaxy is capable of con- 
fining an energetic particle of a given energy and charge, a neces- 
sary precondition in almost all models of particle acceleration. 

The idea that radio-loud AGN might be the origin of the 
UHECR receives some tentative support from, or is at least consis- 
tent with, recent results from the Pierre Auger Observatory (PAO) 
suggesting that the spatial arrival directions of UHECR above 
6 X 10^'' eV are correlated with local AGN (Abraham et al. 2007). 
The imposition of this high low-energy cutoff on the cosmic rays 
ought to imply that they have a relatively local (within ~ 100 
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Mpc) origin, since UHECR at these energies coming from larger 
distances would suffer strong attenuation due to interactions with 
the photons of the cosmic microwave background radiation (the 
so-called Greisen-Zat'sepin-Kuzmin or GZK cutoff; Greisen 1966) 
and also means that these UHECR undergo the smallest possible 
deflection in the Galactic and intergalactic magnetic fields. A par- 
ticularly striking effect in the PAO data released in 2007 was the 
spatial coincidence between several of the UHECR and the position 
of the closest radio galaxy to us, Centaurus A (e.g. Moskalenko et 
al. 2009). While it is not yet clear whether the correlation with local 
AGN remains significant in the PAO data collected since 2007, up- 
dated versions of the Abraham et al. (2007) map appear to show 
a continued overdensity of UHECR around the position of Cen 
A (e.g. Fargion 2009). Meanwhile, several authors have suggested 
that the correlation between the arrival directions of UHECR in the 
original PAO dataset and the positions of local radio-loud AGN is 
at least as good as that with AGN in general (Nagar & Matulich 
2008; Hillas 2009). 

How can specifically radio-loud AGN accelerate UHECR? It 
is of course possible that they are accelerated on sub-parsec scales, 
comparable to the scale of jet generation or initial coUimation. The 
high photon and magnetic field energy densities expected close to 
the active nucleus provide important loss processes, but the acceler- 
ation efficiency might also be higher. Many authors have discussed 
mechanisms by which UHECR can be accelerated in the nuclear 
regions of Cen A and of radio galaxies in general (e.g. KachelrieB, 
Ostapchenko & Tomas 2009) but these necessarily rely on assump- 
tions about the physical conditions close to the nucleus that are 
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hard to test observationally. In what follows I therefore focus on 
the larger-scale components of radio-loud AGN. 

Direct information about the leptonic particle acceleration 
processes in radio galaxies, derived from observations in the op- 
tical and X-ray where the synchrotron loss timescales are shorter 
than the transport timescales from the nuclei so that in situ particle 
acceleration is required, implies that particle acceleration must be 
taking place in the hotspots of powerful double (Fanaroff & Riley 
1974 class II, hereafter FRII) radio galaxies, and in the kpc-scale 
jets of the lower-power FRI class. FRII hotspots have traditionally 
been modelled as the terminal shocks of the relativistic, internally 
supersonic jet that extends up to Mpc scales in these objects (e.g. 
Blandford & Rees 1974; Heavens & Meisenheimer 1987; Meisen- 
heimer et al. 1989), and, while optical and X-ray synchrotron ev- 
idence complicates this picture (e.g. Prieto et al. 2002; Wilson, 
Young Sl Shopbell 2001; Hardcastle et al. 2007a) it seems clear 
that they are particle acceleration sites. Moreover, their sizes and 
their magnetic field strengths, which can be measured very well via 
the inverse-Compton process in the most luminous systems where 
X-ray synchrotron emission is not a contaminant (e.g. Harris et 
al. 1994; Hardcastle et al. 2004) are certainly sufficient to allow 
UHECR to be confined (Hillas 1984). However, the space density 
of FRIIs is very low: we expect only a few within the GZK cut- 
off (for example, the nearest FRII in the northern sky, 3C 98, is at 
a distance of 134 Mpc) and so their effect on the PAO sky above 
6 X 10^^ eV is negligible. 

The numerically dominant population of radio galaxies, by 
several orders of magnitude, within 100 Mpc is composed of low- 
power FRI objects. Here the resolved particle acceleration region is 
typically the 100-pc to kpc-scale inner jet. Several nearby FRI ra- 
dio galaxies, including Cen A (e.g., Hardcastle et al. 2003, 2007c; 
Goodger et al. 2010) and M87 (e.g., Perlman & Wilson 2005; Har- 
ris et al. 2006) have jets that are comparatively strong sources of 
X-ray synchrotron emission, allowing their particle acceleration 
properties to be studied in detail, while the evidence is consistent 
with the idea that all powerful FRI jets can accelerate leptons to 
the > TeV energies required for X-ray synchrotron emission (e.g. 
Worrall et al. 2001). The picture that emerges from the X-ray ob- 
servations is of a combination of strongly localized particle accel- 
eration, which may be due to small-scale shocks, and a more dif- 
fuse process, which produces a different X-ray spectrum (and there- 
fore a different electron energy spectrum) and which may therefore 
have different underlying acceleration physics. It has been argued, 
most recently by Honda (2009), that the Cen A jet is capable of 
accelerating protons to energies comparable to those of the PAO 
UHECR, which of course implies acceleration of heavy nuclei to 
even higher energies. This work relies on rather generous assump- 
tions about the sizes and magnetic field strengths of the acceleration 
regions, though: as yet we have no direct constraint on the magnetic 
field strength in FRI jets (although TeV inverse-Compton emission 
should in principle provide one; Hardcastle &. Croston, in prep.). 

This leaves us with the possibility of UHECR acceleration in 
the lobes, the largest-scale components of both FRI and FRII radio 
galaxies. At first sight these appear less promising candidates for 
UHECR acceleration, since there is little direct evidence for in situ 
particle acceleration in the lobes. However, in the case of the 600- 
kpc giant lobes of Cen A (Hardcastle et al. 2009, hereafter H09) 
we showed that the high-frequency radio data from the Wilkinson 
Microwave Anisotropy Probe {WMAP) are consistent with the idea 
that the lobes contain at least some relatively energetic leptons; 
they do not rule out the idea that particle acceleration is ongoing at 
some level. Similarly, while we do not as yet have a robust inverse- 



Compton measurement of the magnetic field in the lobes of any FRI 
radio galaxy, the available limits in the case of Cen A constrain the 
field strength to be comparable to or greater than the equipartition 
value. H09 argued that the known size, and the limits on B, for 
the giant lobes meant that they could confine protons of energies of 
order 10^" eV, and could therefore accelerate protons to such en- 
ergies, provided that a relatively efficient acceleration process was 
able to operate. We also showed that, provided that the energy in- 
dex for the accelerated cosmic rays is relatively flat, the energetic 
requirements for the acceleration of the PAO UHECR plausibly as- 
sociated with Cen A are trivially satisfied — UHECR need only 
account for a small fraction of the total source energetics. Our pre- 
ferred acceleration mechanism involved scattering off relativistic 
turbulence within the lobes, which requires the assumption that the 
internal energy density is not dominated by thermal particles (see 
also O' Sullivan et al. 2009) but is otherwise consistent with ob- 
servations. We will return to the question of particle content and 
lobe energetics later in the paper, but in the next section I will show 
that a model in which UHECR are accelerated in the giant lobes is 
unique in providing some predictions for the spatial and energetic 
properties of UHECR which may already be testable using the PAO 
data. 

Finally, it should be noted that none of the above mechanisms 
are mutually exclusive. In fact, it seems highly likely that, in a 
source like Cen A, hadronic cosmic rays can be accelerated in the 
nucleus and the kpc-scale jet as well as in the giant lobes. Parti- 
cles accelerated in the inner few kpc will eventually pass into the 
giant lobes and will then be confined (and potentially accelerated) 
there for some time before escaping. Hybrid models of this form 
potentially reduce the problems of acceleration purely in the lobes, 
by providing a seed population of cosmic rays at say 10^^ - lO'^* 
eV and therefore reducing the required UHECR acceleration time 
in the lobes. A corollary of this, unfortunately, is that the ability 
of the giant lobes to confine UHECR, irrespective of whether they 
can accelerate them, implies that the UHECR will be emitted by 
a source like Cen A on scales of the giant lobes, whatever their 
original acceleration site. Even if all UHECR were generated at the 
nucleus, we would not expect a source like Cen A to appear 'point- 
like' at the resolution of the PAO, so we cannot use the observed 
large-scale excess of UHECR around Cen A to argue that accelera- 
tion takes place either wholly or even partly in the giant lobes. This 
limitation should be borne in mind in what follows. 

The remainder of the paper is structured as follows. In Sec- 
tion |2] I show that the requirement that the lobes can confine high- 
energy particles gives a potentially interesting constraint on their 
radio luminosity, and argue that this means that if the PAO UHECR 
are protons they are likely to originate in a small number of bright 
nearby radio galaxies, all probably nearby well-studied objects. In 
Section[3]l discuss our best existing constraints on the particle con- 
tent of FRI lobes and the implications for cosmic ray acceleration 
and composition. Finally, in Section|4]l discuss the implications of 
a picture in which particles are accelerated in radio galaxy lobes 
for the future of cosmic ray astronomy. Throughout the paper I use 
a cosmology in which Hq = 70 km s^^ Mpc^^, Sim = 0.3 and 
SIa = 0.7. The distance to Cen A is taken to be 3.7 Mpc (the mean 
of 5 distance estimates given in Ferrarese et al. 2007). 

2 CONSTRAINTS ON THE RADIO EMISSION 

In H09 we argued that stochastic acceleration by large-scale mag- 
netic turbulence (as discussed by, e.g., Stawarz & Petrosian 2008) 
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imposes a condition on the energy, acceleration region radius and 
magnetic field strength that is equivalent to the classical particle 
confinement condition. This comes about because the acceleration 
timescale, assuming Bohm diffusion, is ~ lOri^/c, where tl is 
the Larmor radius, while the timescale for diffusive escape from 
the lobes is ~ 3R^/ri^c: equating these two gives ri, ^ R for 
the highest-energy cosmic rays that can be accelerated efficiently, 
which is simply the confinement condition. Ignoring numerical fac- 
tors of order unity, therefore, we can consider the confinement con- 
dition as giving us the (best-case) estimate of the upper limit on the 
cosmic ray energy. In this section I demonstrate that the require- 
ment that the lobes be capable of confining UHECR at the energies 
observed gives rise to an interesting constraint on the combination 
of the radio luminosity and size of the lobes. 

The confinement criterion for a particle of energy Ep is that 
the gyroradius tl be less than the size of the region R; in other 
words, in SI units. 
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where Z is the nuclear charge and e is the charge on the proton. 
We consider a spherical radio lobe with radius R and a uniform 
magnetic field strength and electron energy density. Let the electron 
energy distribution be given by N{Ee) and let the magnetic field 
be a factor e away from equipartition, so that 
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Here, as in H09, we are assuming 'true' equipartition between the 
electron energy spectrum and the magnetic field, by integrating 
over all electron energies [following Myers & Spangler (1985) and 
Hardcastle, Birkinshaw & Worrall (1998)] rather than between the 
energies corresponding to a pair of observed frequencies as in 'clas- 
sical' equipartition. The differences between the two equipartition 
formulae are discussed in more detail by Brunetti, Setti & Comas- 
tri (1997) and Beck & Krause (2005). Beck & Krause (2005) show 
that the 'classical' formula can lead to a significant underestimate 
of the field strength, and thus to the ability to confine high-energy 
particles, in radio galaxy lobes. 

The simplest electron energy distribution we can consider is 
then a power law in energy, i.e. N{E^) = NqE'", with a mini- 
mum and maximum energy given by iSmin and iSmax respectively. 
This allows us to solve the integral of eq.|2]analytically. Let 
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(3) 

The total energy in electrons is then A'^o^- However, shock- 
acceleration models predict that the electron energy spectrum 
should actually be described by a power law in momentum, i.e. 
n{p) = nop^" (e.g. Blandford & Ostriker 1978). In this case, the 
electron energy integral becomes 

1/2 



1 + 



p-'dp (4) 



(where Pmin and pmax are the momenta corresponding to the ener- 
gies i?min and i?max respectively, and the leading factor accounts 
for the difference between the normalizations A'^o and no in energy 
and momentum). It is most convenient to evaluate the integral of 



equation |4] numerically, though clearly it converges to the analyt- 
ical solutions of equation [3] in the limit that i?min 2> me(? . We 
comment below on the differences that arise when using these two 
values of /. 

Now for a power-law electron energy distribution the vol- 
ume emissivity in synchrotron emission at a given (rest-frame) fre- 
quency V may be written 



J{u) = C{s)Nou 
(Longair 1994 eq. 18.49) where 
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for an isotropic pitch angle distribution. The frequency dependence 
in eq. |5] expresses the well-known relationship between s and the 
synchrotron spectral index ct. This result is valid both for the trun- 
cated power-law distribution used in equation [3] and for the power 
law in momentum described in equation |4] so long as the chosen 
observing frequency lies in a region where the electron spectrum 
is a power law (i.e. not too close to the high-energy or low-energy 
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we can now use equation |7] to eliminate the electron spectral nor- 
malization A'o from equation[5] 
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Now rewriting eq. [T]as B > Ep/ZeRc, we can eliminate B from 
equation[8l turning it into an inequality: 
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Finally, we can turn the emissivity J{v) into an observable quantity 
by noting that L{iy) = |7ri?^ J(i^), so that 
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We have derived a limit on the luminosity and size of a lobe that 
is (marginally) capable of confining a particle of energy Ep and 
charge Z. For conventional values of s (in the range 2-3) note the 
strong dependence of the luminosity on Ep/Z (rigidity), the linear 
dependence on the equipartition factor e, and the relatively weak 
dependence on source size R, which is in the sense that a lower lu- 
minosity is required for a larger size. For known e, s and Ep/Z, eq. 
|9] defines a line in the conventional radio luminosity/size diagram 
for radio galaxies separating those that can accelerate such particles 
from those that cannot (Fig.[T), if we adopt a model such as that of 
H09 in which efficient stochastic particle acceleration is possible in 
the lobes. 

What are the implications for the population of radio galaxies 
that can accelerate UHECR? First of all, we can substitute physical 
constants into eq. |9]to obtain a relationship in useful units. If we 
take s = 2, E^in = 5 MeV, E'max = 5 GeV, u = 408 MHz, 
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Figure 1. The luminosity-size constraint of eq. [TO]for e = 1, E20/Z = 
1 applied to both lobes (i.e. D = 4iJ, = 2L) and plotted 

over the luminosity-size diagram for the B2 bright sample (LAS data 
from Fanti et al. 1987, flux densities, redshifts and spectral indices 
as tabulated by Hardcastle et al. 2003) and 3CRR sample (data from 
|http ://3crr.extragalactic.info/ ). Only sources above the solid line can con- 
fine UHECR with energies of 10^'' eV; sources below the line cannot. 
While almost all the powerful 3CRR sources can, in principle, many of 
the lower-luminosity B2 sources cannot. Note that this figure is illustrative 
only; only a few of the lowest-luminosity objects in this plot are within the 
GZK cutoff and none are in the southern sky to which the PAO is most sen- 
sitive. No attempt has been made to take projection into account or to model 
the actual physical sizes of the lobes. 



R ~ Tioo X 100 kpc, E — E20 X 10^" eV, and use the numerically 
calculated expression for I based on a power-law distribution in 
momentum, then we obtain 
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What restrictions does this put on the population of radio galaxies 
capable of accelerating UHECR? We can begin by turning this into 
a strict limit on luminosity by imposing the observationally-based 
limit that R < 250 kpc (i.e. rioo < 2.5) since we know that very 
few radio galaxies exceed 1 Mpc in size. To compare to total radio 
luminosity we must also scale up by a factor 2, since so far we 
have only been considering the luminosity of a single lobe. This 
gives us a strict lower limit on I/408 of 2.5 x 10^* W Hz~^ for 
e — 1.0, E20 = 1, Z — 1. Immediately we see (Fig.lTJ that only 
reasonably luminous radio galaxies can accelerate UHECR to these 
energies; the Fanaroff-Riley break is at ~ 3 x 10^^ W Hz"^ at 408 
MHz. (Centaurus A, with a 408-MHz luminosity ~ 3 x 10^" W 
Hz~^, just satisfies this criterion, as we would expect given the 
results of H09.) If we compare with a recent determination of the 
local radio luminosity function, such as that by Mauch & Sadler 
(2007), we find that within a sphere of radius 100 Mpc (the right 
order of magnitude for 10^" eV protons) we expect ^ 20 radio 
galaxies satisfying the luminosity criterion alone, of which not all 
will satisfy the size criterion (and of course only roughly half of 
which will be visible to the PAO). Thus for the parameters we have 
used here we see that radio galaxies capable of accelerating protons 
to 10^° eV will be rare. In addition, since their luminosities are 
large and their distances constrained, their fluxes are known (> 2 
Jy at 408 MHz) and so we can say that all such objects are probably 
bright, well-studied local radio galaxies. 




Figure 2. The normalizing luminosity of equation [To] as a function of the 
power-law index s of the energy/momentum power law. The lower (black) 



curve shows E„ 



0.5 MeV. the middle (red) curve shows Ej^ 



MeV. and the upper (green) curve shows -Emin = 50 MeV. The solid lines 
show the luminosity for an assumed power-law distribution in momentum 
and the dashed lines show a power law in energy. The clearest trend is a 
decrease in the limiting luminosity with increasing s. We also see that the 
effect of changing E^i^ is very limited for s = 2 but very significant by 
the time s = 3. The effect of incorrectly assuming a power-law distribution 
in energy is only significant for the lowest value of -Emin- 



Let us now consider varying some of the assumptions in the 
calculations above. 



2.1 Power-law index and minimum energy 

The normalizing luminosity in equation [TO] has a relatively strong 
dependence on the power-law index s. This reflects the fact that 
the electron energy density, and thus the magnetic field strength, 
is dominated by the low-energy electrons, while it is high-energy 
electrons that produce the observed radio emission at our normaliz- 
ing frequency. There is also a dependence on the minimum energy 
Emin which is stronger for larger s. These dependencies are illus- 
trated in Fig. 12] Values of s close to 2.0 are predicted in shock ac- 
celeration models and appear to be consistent with observation (cf. 
Young et al. 2005). Fig.|2]shows that the normalizing luminosity is 
not greatly affected for values close to 2, say s = 2.2, so that the 
number of potential accelerating sources is probably not greatly af- 
fected by our uncertainties on this parameter or on the appropriate 

value of Emin- 

The dependence on iSmax is always weak, and so it is unnec- 
essary to put in a more realistic electron energy distribution with 
spectral steepening below i?max- 



2.2 Particle energy 

The number of radio galaxies capable of accelerating UHECR pro- 
tons to high energies is a very strong function of £20- Even for 
E20 = 2, the expected number of radio galaxies in the southern 
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sky within 100 Mpc that satisfy eq. [TO] neglecting the size con- 
straint, is less than 1 (and here we also neglect the steep decrease 
in the appropriate radius to use due to the energy-dependent GZK 
cutoff). By contrast, if we set E20 ~ 0.6, there are perhaps 40 radio 
galaxies in the southern sky that are in principle capable of acceler- 
ating protons to those energies, again neglecting the size constraint. 
Effectively, therefore, this model for UHECR acceleration predicts 
a very steep cutoff in the integrated source spectrum of UHECR 
protons which, by chance, occurs at energies close to the energy at 
which GZK effects become significant, and which therefore rein- 
forces the effect of the GZK cutoff. 



2.3 Equipartition 

In the calculations above I have used e — 1, corresponding to 
equipartition between radiating particles and magnetic field, which 
is consistent with the known constraints from inverse-Compton ra- 
diation from Cen A. If there were an energetically dominant popu- 
lation of non-radiating charged particles, such as protons, we would 
expect the magnetic field to be in equipartition with those and so 
e ^ 1. On the other hand, the evidence in those FRII sources in 
which inverse-Compton modelling has been possible is that the 
magnetic field strength is typically somewhat below the equipar- 
tition value, implying e > 1. Values of e 2> 1, implying very low 
B-field strengths for a given observed synchrotron luminosity, ob- 
viously make it very hard for lobes to accelerate UHECR. Values 
<^ 1 make it easier, but given the rather flat radio luminosity func- 
tion do not immediately fill the sky with UHECR-emitting radio 
galaxies: for example, substituting e = 0.1, E20 = 1, Z = 1, 
noo = 2.5 in eq.[TO]and integrating over the luminosity function 
gives around 50 radio galaxies in the southern sky that meet the 
luminosity constraint. 



2.4 Relativistic turbulence 

The acceleration timescale and therefore the efficiency of UHECR 
acceleration in the acceleration model of H09 depend on the 
presence of strong magnetic turbulence in the lobes: we require 
CA;urb ~ Uo where Utmh and Uo are the energy densities in the 
turbulent magnetic field component and the unperturbed compo- 
nent respectively. Once (Aurb becomes much less than Uo turbu- 
lent acceleration will be very much less efficient. Magnetic turbu- 
lence is presumably generated by large-scale hydrodynamic pro- 
cesses and therefore relies on a continued energy supply by the jet. 
In lobes where the jet is disconnected the turbulence will decay on 
timescales which may be as short as R/c. Since the energy density 
in relativistic turbulence may therefore depend on local details of 
the coupling between the jets and the large-scale lobes, it is clearly 
therefore possible to imagine a situation in which the particle accel- 
eration efficiency varies from source to source, depending on such 
factors as the large-scale morphology of the lobes and the presence 
or absence of strong jet-lobe interactions (such as jet termination 
shocks). We necessarily cannot take account of this in the simple 
models presented here, but it should be borne in mind that the lu- 
minosity/size cutoffs apply only in sources in which relativistic tur- 
bulence can efficiently be maintained; in particular, disconnected 
lobes (possibly even including the S giant lobe of Cen A; see H09) 
are likely to be unable to accelerate UHECR. 



2.5 Composition 

As with all UHECR acceleration models, it is much easier to ac- 
celerate nuclei, with Z > 1, than protons to a given energy. If we 
naively substitute e = 1.0, -E20 = 1, Z — 26 (iron), rioo = 2.5 
into eq. (TO] then the luminosity limit comes down by nearly 5 or- 
ders of magnitude and we find that practically every radio galaxy 
in the sky is a potential UHECR source, although in practice the 
size constraint will still impose some limitations. A self-consistent 
model for heavy nucleus acceleration in the lobes would need to 
take account of losses to photodisintegration within the lobes them- 
selves — both the acceleration and loss timescales are shorter for 
nuclei than for protons — but such a model is beyond the scope 
of this paper. Here we have simply to note that the strong rigid- 
ity dependence of eq.[9]means that the composition of the baryons 
available to accelerate in the source will have a strong effect on 
the predicted composition, energy spectrum and arrival positions 
of UHECR in this model (all of course modified by propagation 
effects; e.g. Hooper & Taylor 2009). We discuss the available con- 
straints on composition in the next section. 



3 PARTICLE CONTENT 

The preceding sections have shown a strong dependence of the 
predictions of a model in which UHECR are accelerated in radio 
galaxy lobes on the source composition of the particles in the lobes, 
and so at this point it is appropriate to comment on the known con- 
straints on the particle content of lobes, and to ask what sources of 
(1) protons and (2) heavy nuclei are available on these scales. 

We do not know whether jets in radio-loud AGN are electron- 
proton or electron-positron in their initial composition. There is 
some evidence in FRII radio galaxies that the lobes are not dom- 
inated energetically by protons — see Croston et al. (2005) for the 
argument — but there is no way in these FRII lobes to rule out the 
possibility that a relativistic proton population has energy roughly 
comparable to that in the electrons and magnetic field, whose en- 
ergy densities can be measured. A fortiori, we do not know the 
expected fraction of heavy nuclei in these lobes, since even if the 
jets contain protons, it is not clear at what point in the jet generation 
process they get there. 

The situation in FRI radio galaxies has been known for many 
years to be more interesting. Here inverse-Compton measurements 
are not in general available, but the minimum pressures (approxi- 
mately equivalent to the assumption of equipartition between field 
and electrons alone) in the large-scale lobes or plumes can be sev- 
eral orders of magnitude below those of the external medium (see 
Hardcastle et al. 2007b and references therein) and, since inverse- 
Compton constraints rule out electron dominance by very large fac- 
tors and energetic dominance by magnetic field seems a priori im- 
plausible, it is conventional to suggest that the missing pressure is 
supplied by a population of non-radiating particles (perhaps with 
the magnetic field in approximate equipartition, implying e <C 1). 
We have recently argued (Croston et al. 2008; Croston & Hardcas- 
tle, in prep.) that there is some evidence that these particles are the 
same particle population that is required to be entrained to decel- 
erate the kpc-scale jet; this requires a means of efficiently heating 
or accelerating these particles to make them provide the required 
pressure. Entrainment is interesting here because it does give us 
some constraints on the expected abundance of the heavy particle 
population, which should be similar to that of the external medium 
(i.e. roughly 1 iron nucleus per lO"" protons for an assumed 0.3 
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solar abundance). It is important to note, though, that if the en- 
trained particles were thermal and dominated the energetics, then 
the Alfven speeds in the lobes would be ^ c and the efficiency of 
stochastic acceleration in the lobes would be greatly reduced (H09; 
O'Sullivan et al. 2009). To avoid this, we would need the entrained 
particles to be relativistic and to participate in the equipartition pro- 
cess so that the energy densities in magnetic field and baryons were 
comparable. There is as yet no direct evidence that rules this pic- 
ture out (see Croston & Hardcastle, in prep., for more discussion of 
constraints on the state of the entrained material). In addition, we 
have shown that the necessity for entrainment varies from source to 
source, even among FRIs (Croston et al. 2008). It is certainly possi- 
ble that jets in the low-power sources are initially electron-positron 
and that the amount of material entrained, at least in some sources, 
is enough to provide the seed population for stochastic accelera- 
tion of baryons in the lobes while not being so much that stochastic 
acceleration is inefficient. 



4 PREDICTIONS AND THE OUTLOOK FOR COSMIC 
RAY ASTRONOMY 

One of the attractive features of the model proposed in H09 and 
discussed in this paper is that, at least superfically, it makes some 
simple testable predictions. If the UHECR mapped by the PAG are 
protons, which would be implied by a genuine detection of cor- 
relation on the sky with the positions of distant objects (perhaps 
excluding Cen A; see below) then I have shown above that the lobe 
acceleration model would imply that their local sources are phys- 
ically large (> 100 kpc), luminous, relatively rare radio galaxies. 
A couple of hints that this is so are already seen, firstly in the ap- 
parent excess of events around Cen A, and secondly in the work 
of Nagar & Matulich (2008) who found a correlation between the 
positions of extended radio galaxies and the arrival positions of the 
R^O UHECR. With the eventual release of updated PAG positional 
data it should be possible to make a systematic investigation of all 
possible radio-loud sources of UHECR. 

However, the situation is complicated by the current uncer- 
tainty about the composition of the PAO cosmic rays. Current mea- 
surements of the mean and RMS depth of shower maximum imply 
a large fraction of heavy or intermediate mass nuclei, and certainly 
do not appear consistent with a pure-proton spectrum (Abraham et 
al. 2010), although it is important to note that the HiRes results 
are quite different (e.g. Aloisio, Berezinsky & Gazizov 2009). If 
we are required to accelerate nuclei with Z > 1, then the predic- 
tions change in two crucial ways: firstly, the number of potential 
radio-galaxy UHECR sources increases rapidly with increasing Z, 
as discussed above (Section |231 l; secondly, it becomes increasingly 
unlikely that a spatial correlation will be observed between the ar- 
rival directions of the UHECR and the positions of their sources, 
due to the larger deflection of Z > 1 UHECR in intergalactic mag- 
netic fields. At this point a model that intended to reproduce the 
observations would need to take account of (1) the spatial distribu- 
tion of the radio galaxy sources throughout the GZK volume and 
perhaps beyond; (2) the distribution of the physical conditions in 
their lobes; (3) the intrinsic UHECR energy and composition spec- 
trum; (4) propagation losses for the various species of UHECR; and 
(5) deflection in the Galactic and intergalactic magnetic fields. The 
tools to do (I) and (2) are available and to some extent presented 
in this paper; (3) remains very uncertain, though we have some 
constraints (see Section |3] above); and (4) and (5) are in principle 
possible (e.g.. Hooper & Taylor 2009), although crucial elements 



remain uncertain. Putting all five elements together would be a ma- 
jor effort which may not yet be justified by the state of the data, 
but detailed modelling like this will be the way forward if we are to 
start doing serious astrophysics with UHECR observations. 

Having said this, some limiting cases of the model if the PAG 
UHECR are heavy or intermediate-mass nuclei are relatively easy 
to imagine. Aloisio et al. (2009) discuss what they call the 'disap- 
pointing model' for the PAG results in which rigidity-dependent ac- 
celeration (as is implicit in the lobe-acceleration picture) together 
with an acceleration cutoff for protons at ~ 10^^ eV leads to a 
steadily increasing fraction of heavy nuclei with increasing energy 
across the PAG band, as observed. Given that we require the most 
favourable assumptions to make radio galaxies accelerate protons 
up to 10^" eV, it is easy to imagine that the high-energy cutoff might 
be reduced by an order of magnitude or so (e.g., by reducing the 
maximum energy placed into large-scale magnetic turbulence), so 
that they might provide the 'disappointing' population required by 
Aloisio et al.The detection of an enhanced count rate around Cen 
A would then be explained by its proximity, which has the effect 
that few-Z particles are deflected by the Galactic and intergalac- 
tic magnetic fields by only a few (up to 10) degrees. However, the 
UHECR from all other sources would be scattered by much larger 
angles and it would never be possible to identify them with their 
parent radio galaxies. While this model would be slightly less dis- 
appointing than the limiting case suggested by Aloisio et al., Cen 
A would remain the only detectable UHECR source in the sky. 



5 SUMMARY 

The principal results of this paper may be summarized as follows: 

(i) Stochastic acceleration of UHECR in the large-scale lobes of 
radio galaxies may be possible, but there are strong (though model- 
dependent) constraints on the properties of the radio galaxies that 
can accelerate them to the highest energies (10^" eV). 

(ii) These constraints imply that only a small number of local 
radio galaxies can be involved in the acceleration of UHECR, if 
the UHECR are protons, and that UHECR energies will cut off 
steeply around the energies currently being observed by the PAG; 
this model is testable in principle using existing radio surveys and 
up-to-date UHECR arrival positions, and is consistent with much 
of the available data. 

(iii) However, if UHECR are heavy nuclei with Z > 1, as sug- 
gested by the latest PAG composition results, then many more radio 
galaxies can be sites of UHECR acceleration, and it may be that the 
nearest radio galaxy, Cen A, will be the only identifiable source in 
the cosmic-ray sky. 



ACKNOWLEDGEMENTS 

This work was first presented at the Trondheim workshop on 
'Searching for the origins of cosmic rays' in June 2009, and I am 
indebted to the organizers for inviting me and to many participants 
there for helpful comments. I also owe a debt of gratitude to Lukasz 
Stawarz and Teddy Cheung, without whom the discussion of cos- 
mic ray acceleration in H09 would have been extremely limited. 
The paper was substantially improved as a result of comments from 
an anonymous referee. I acknowledge generous financial support 
from the Royal Society through the University Research Fellow- 
ships scheme. 



Radio galaxies and cosmic rays 1 

REFERENCES 

Abraham, J., et al. [for the Pierre Auger Collaboration], 2007, Sci, 318, 938 
Abraham, J., et al. (for the Pierre Auger Collaboration), 2010, Phys. Lett. 

B., in press (arXiv: 1002. 1975) 
Aloisio, R., Berezinsky, V., Gazizov, A., 2009, arXiv:0907.5194 
Beck, R., Krause, M., 2005, Astron. Nachr., 326, 414 
Blandford, R.D., Ostriker, J.R, 1978, ApJ, 221, L29 
Blandford, R.D., Rees, M.J., 1974, MNRAS, 169, 395 
Brunetti, G., Setti, G., Comastri, A., 1997, A&A, 325, 898 
Burbidge, G., 1956, ApJ, 124, 416 

Croston, J.H., Hardcastle, M.J., Harris, D.E., Belsole, E., Birkinshaw, M., 

Worrall, D.M., 2005, ApJ, 626, 733 
Croston, J.H., Hardcastle, M.J., Birkinshaw, M., Worrall, D.M., Laing, 

R.A., 2008, MNRAS, 386, 1709 
Fanaroff, B.L., Riley, J.M., 1974, MNRAS, 167, 31P 
Fanti, C, Fanti, R., De Ruiter, H.R., Parma, P, 1987, A&AS, 69, 57 
Fargion, D., 2009, arXiv:0908.2650 

Ferrarese, L., Mould, J.R., Stetson, P.B., Tonry, J.L., Blakeslee, J.P., Ajhar, 

E.A., 2007, ApJ, 654, 186 
Goodger, J.L., et al., 2010, ApJ, 708, 675 
Greisen, K., 1966, Phys. Rev. Lett., 16, 748 

Hardcastle, M.J., Birkinshaw, M., Worrall, D.M., 1998, MNRAS, 294, 615 

Hardcastle, M.J., Birkinshaw, M., Cameron, R., Harris, D.E., Looney, L.W., 
Worrall, D.M., 2002, ApJ, 581, 948 

Hardcastle, M.J., Worrall, D.M., Kraft, R.P., Forman, W.R., Jones, C, Mur- 
ray, S.S., 2003, ApJ, 593, 169 

Hardcastle, M.J., Harris, D.E., Worrall, D.M., Birkinshaw, M., 2004, ApJ, 
612, 729 

Hardcastle, M.J., Croston, J.H., Kraft, R.P, 2007a, ApJ, 669, 893 
Hardcastle, M.J., Kraft, R.P, Worrall, D.M., Croston, J.H., Evans, D.A., 

Birkinshaw, M., Murray, S.S., 2007b, ApJ, 662, 166 
Hardcastle, M.J., et al., 2007c, ApJ, 670, L81 

Hardcastle, M.J., Worrall, D.M., Birkinshaw, M., Canosa, CM., 2003, MN- 
RAS, 338, 176 

Hardcastle, M.J., Cheung, C.C., Feain, I.J., Stawarz, L., 2009, MNRAS, 
393, 1041 

Harris, D.E., Carilli, C.L., Perley, R.A., 1994, Nat, 367, 713 

Harris, D.E., Cheung, C.C., Biretta, J.A., Sparks, W.B., Junor, W., Perlman, 

E.S., Wilson, A.S., 2006, ApJ, 640, 211 
Heavens, A.F., Meisenheimer, K., 1987, MNRAS, 225, 335 
Hillas, A.M., 1984, ARA&A, 22, 425 
Hillas, A.M., 2009, Astropart. Phys., 32, 160 
Honda, M., 2009, ApJ, 706, 1517 
Hooper, D., Taylor, A.M., 2009, Astropart. Phys., 33, 151 
KachelrieB, M., Ostapchenko, S., Tomas, R., 2009, NJPh 11 065017 
Longair, M.S., 1994, High energy astrophysics, Cambridge University 

Press, Cambridge 
Mauch, T., Sadler, E., 2007, MNRAS, 375, 931 

Meisenheimer, K., Roser, H.-J., Hiltner, PR., Yates, M.G., Longair, M.S., 

Chini, R., Perley, R.A., 1989, A&A, 219, 63 
Moskalenko, I.V., Stawarz, L., Porter, T.A., Cheung, C.C., 2009, ApJ, 693, 

1261 

Myers, S.T., Spangler, S.R., 1985, ApJ, 291, 52 
Nagar, N.M., MatuHch, J., 2008, A&A, 488, 879 
O'SuUivan, S., Reville, B., Taylor, A.M., 2009, MNRAS, 400, 248 
Perlman, E.S., Wilson, A.S., 2005, ApJ, 627, 140 
Prieto, M.A., Brunetti, G., Mack, K.H., 2002, Sci, 298, 193 
Stawarz, L., Petrosian, V., 2008, ApJ, 681, 1825 
Wilson, A.S., Young, A.J., Shophell, PL., 2001, ApJ, 547, 740 
Worrall, D.M., Birkinshaw, M., Hardcastle, M.J., 2001, MNRAS, 326, L7 
Young, A., Rudnick, L., Katz, D., DeLaney, T, Kassim, N.E., Makishima, 
K., 2005, ApJ, 626, 748 



